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The Floating Potential Measurement Unit (FPMU) has detected high negative amplitude 
rapid charging events (RCEs) on the International Space Station (ISS) at the morning 
terminator. These events are larger and more rapid than the ISS morning charging events 
first seen by the Floating Potential Probe (FPP) on ISS in 2001. In this paper, we describe a 
theory for the RCEs that further elucidates the nature of spacecraft charging in low Earth 
orbit (LEO) in a non-equilibrium situation. The model accounts for all essential aspects of 
the newly discovered phenomenon, and is amenable to testing on-orbit. Predictions of the 
model for the amplitude of the ISS RCEs for the full set of ISS solar arrays and for the 
coming solar cycle are given, and the results of modeling by the Environments WorkBench 
(EWB) are compared to the observed events to show that the phenomenon can be explained 
by solar array driven charging. The situation is unique because the coverglasses have not yet 
reached equilibrium with the surrounding plasma during the RCEs. Finally, a prescription 
for further use of the ISS for investigating fundamental plasma physics in LEO is given. 
Already, plasma and charging monitoring instruments on ISS have taught us much about 
spacecraft interactions with the dense LEO plasma, and we expect they will continue to yield 
more valuable science when the Japanese Experiment Module (JEM) is in place. 
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K 

= 

dielectric constant of a material 


= 

Debye length (m) 

m e 

- 

electron mass 

N e 

= 

electron density (m' 3 ) 

NLP 

= 

Narrow (Sweep) Langmuir Probe 

PCU 

= 

Plasma Contacting Unit 

0 

= 

normal charging peak value (Volts) 

<P 

= 

theoretical RCE amplitude 

<Pflt 

= 

RCE amplitude from curve fits 

q or Q 

= 

charge (Coulombs) 

RCE 

= 

Rapid Charging Event 

r sh 

= 

sheath radius (m) 

T 

1 e 

= 

electron temperature (eV or Kelvins) 

V 

= 

ISS velocity 

vxB’l 

= 

induced potential from magnetic field 

V 

= 

potential with respect to the plasma (V) 

Vobs 

= 

observed ISS structure charging (Volts) 

V pk 

= 

peak voltage (V olts) 

WLP 

= 

Wide (Sweep) Langmuir Probe 


I. History of ISS Charging Theory 

T HE original expectation (1991). In 1990 and 1991, a Tiger Team was convened by the Space Station Freedom 
(SSF) program to determine whether SSF would charge up due to its high voltage solar arrays and its negatively 
grounded power system. The report of the Tiger Team was that SSF would likely charge up to about -140 V with 
respect to the plasma, that this voltage would break down its dielectric thermal control coatings, and that the 
potential should be actively controlled l ’ 2 . This led to the development and adoption of the PCUs (plasma contactor 
units) that have now been in place on ISS from the year 2000, and are still routinely used to prevent charging when 
EVAs are in progress. 

A. Revision due to electron choking (1991-1994) 

During the Tiger Team activities, it was discovered both analytically 3,4 and experimentally (through ground tests 
in plasma chambers 5 ) that electron collection on the ISS solar cells is effectively choked offby the potentials on the 
coverglass surfaces and in the gaps between solar cells for electron temperatures above about 0.5 eV, and for fully 
neutralized coverglass surfaces. For electron temperatures less than that, electron collection (and the amount of ISS 
charging) would increase strongly with decreasing electron temperature. 

1. SAMPIE vs. PASP+ 

Two flight experiments were performed in 1993 and 1994. The Solar Array Module Plasma Interactions 
Experiment (SAMPIE), a flight experiment conducted on the Space Shuttle, used ISS-type solar cells cemented to a 
solid substrate, and seemed to show that electron choking would be very effective at electron temperatures of about 
0.2 eV 6 . The other experiment, Photovoltaic Array Space Power - Plus (PASP+), performed on a Pegasus-launched 
satellite with ISS-type cells on a spring-loaded flexible substrate, seemed to show that electron choking was not 
nearly as effective as SAMPIE had implied 7 . Ground tests 8 were interpreted as indicating that the flexible substrate 
tests would be more representative of collection on the real ISS arrays, so models were adjusted to simulate PASP+ 
like collection for ISS arrays. After the FPP was installed on ISS in 2000, models were re-adjusted back to agree 
with the SAMPIE results, as no agreement with the measured charging could be obtained with the PASP+ values. 

B. Revision due to extra ion collection (2002, 2006) 

Post-2000 models also could only be made to fit the FPP observations if two more sources of ion collection were 
included. These were ion collection from the (grounded) solar array mast wires and another mysterious conducting 
area probably located on the Russian segment 34,35,36 . 
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1. Russian blanket testing 


In order to determine the ion collection on the Russian segments, samples of Russian thermal blankets were 
tested in plasma chambers at the Marshall Space Flight Center 9 . It was found that due to conducting material 
incorporated in the blankets, some 20-35 square meters of ion collecting area existed on the Russian segment, and 
this would also lead to smaller than expected charging. Incorporating both of these ion collection areas led to very 
good agreement of physical models with the FPP-observed charging values 36 . 

C. Application to “normal” charging 

A steady-state physical model for LEO charging called the Environments WorkBench (EWB) model 10,37 was 
developed by Science Applications International Corporation (SAIC). A specialized version of EWB called the 
Plasma Interactions Model (PIM) was jointly developed by SAIC and Boeing that agreed in all essential aspects 
with the FPP observations. The original PIM code included electron current collection models for the solar arrays 16 
as well as collection on solar arrays, masts, and vehicle structure with a distribution of potential due to vxB»l 

34,35,36,38 (see below ) 


II. Previous Observations of Charging 

A. The original FPP - capabilities and limitations 

The original Floating Potential Probe (FPP) was placed into service on ISS in December, 2000. It had two 
probes, a Langmuir probe to determine plasma characteristics and the structure potential if less than about 1 0 volts, 
and a floating potential probe to measure structure potential between 0 and -150 V. For structure potentials of 
magnitude greater than about 10 volts, no plasma data could be obtained. Although data were taken on floating 
potential every 0.1 second, they were analyzed in 20 second increments (the duration of a single Langmuir probe 
sweep). FPP operated for a period of about 4 months until April, 2001. For a full account of the original FPP, see 
Ferguson". 

B. “Normal” charging amounts and characteristics 

1. Charging 

Charging as seen by FPP occurred mainly near eclipse exits, when the active side of the array was ram-facing, 
before most solar array strings had been turned off for power management, and when the plasma temperature was 
still low. Ferguson and Morton 12 showed that the normal charging lasted for several minutes after eclipse exit and 
could be described in terms of plasma parameters by the equation: 

O = -2.6907 N e 0 1057 e' 8 02Te , (1) 

where N e is the electron density and T e is the electron temperature at eclipse exit. The steep dependence of this 
charging on electron temperature is reflective of the electron-choking effect mentioned earlier. 

For the single solar array in place when FPP was operational, no charging of magnitude above about 30 volts 
was observed, and the charging peaks typically lasted for a few minutes after eclipse exit. Ferguson and Gardner 13,14 
extrapolated the charging to a fUlly array-populated ISS and found that normal charging in excess of 80 volts could 
occur. 

2. vxB*l 

Tthe electric field induced by the ISS motion through the ambient magnetic field, the vxB*l effect, can be 
important. At its maximum value, it can produce a difference of about 0.3 volts/m in potential along a truss 
member. For a 100 meter long truss on ISS, vxB«l can induce up to a 30 V difference in local potentials relative to 
the plasma potential 15 . This effect was quite well modeled by PIM (see for example Mandell et al. 16 ), and can lead 
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to locally heightened or lessened ion or electron current collection. In what follows, we will take out this effect 
when necessary, and will deal only with potentials due to the difference in electron and ion collection of the ISS 
solar arrays and structure. 


m. FPMU 

A. Description 

The FPMU is a package of four plasma instruments designed and built by Utah State University (USU) under 
contract to the NASA/Johnson Space Center (JSC). The purpose of the FPMU is to measure the local potential of 
the ISS relative to the plasma and to measure the density and electron temperature of the local plasma. The four 
probes are: a Floating Potential Probe (FPP), a Plasma Impedance Probe (PIP), a Wide-sweep Langmuir Probe 
(WLP), and a Narrow-sweep Langmuir Probe (NLP) with associated electronics. The operation is autonomous, with 
either an on or an off state. The FPMU has been described elsewhere [Swenson et al. 27,28,29 , Wright et al. 25 ] but a 
brief description of each probe is given below. Initial FPMU on-orbit results are given in Wright et al. 25 , Baijatya et 
al. 33 , and Coffey et al. 26 . 

1. Floating Potential Probe (FPP) 

The FPP is a gold-plated sphere of radius 5.08 cm and is essentially a high impedance (~ 10 11 Q) voltmeter with 
the other connection referenced to ISS ground. The sphere potential is determined by local plasma conditions, and 
within a few kT e of the plasma potential, provides a reference for measuring the potential of the ISS (k = 
Boltzmann’s constant). Data is sampled at 128 Hz with 100 mV resolution. This is the primary probe for ISS 
floating potential because it is a direct measurement of the floating potential. 

2. Wide-Sweep Langmuir Probe (WLP) 

The WLP is a gold-plated sphere also of radius 5.08 cm. The voltage applied to the probe is swept from -20 V to 
+80 V relative to chassis ground (ISS structure) and the resulting currents to the probe are measured. Sweeps are 
accomplished each second, with the potential sweeping from low to high voltage in one second and back down from 
high to low in the next second. The sweep is comprised of three parts as described by Wright et al. 25 [2008] and 
references therein. Step sizes of ~ 25 mV from 0 V to 50 V provide sufficient resolution for a determination of T e 
(which requires several samples in the electron retarding portion of the sweep). The floating potential can be 
obtained over the full -20 V to +80 V range, within an uncertainty of +/- 2 V. The WLP is the primary probe for 
densities because its wide voltage sweep allows complete coverage of the ion and electron saturation regions over a 
wide range in floating potential. Plasma temperatures derived from the WLP are considered secondary to those 
derived from the NLP described below. 

3. Narrow-Sweep Langmuir Probe (NLP) 

The NLP is a guarded gold-plated cylinder with collector dimension of radius 1.43 cm and length 5.08 cm. The 
NLP is placed mid-way on the boom supporting the FPP. A sweep from - 4.9 V to + 4.9 V, in steps of ~12 mV, is 
applied to the NLP during one second, followed by a sweep down from 4.9 V to - 4.9 V in the next second. This 
sweep voltage is referenced to the floating potential measured by the FPP. In this way, the electron and ion 
saturation regions of the plasma current-voltage profile are seen, as the region sampled will move through the 
- 180 V to + 180 V range of the FPP. This configuration allows the density and electron temperature to be 
determined at 1 Hz. Electron temperature derived from this probe is considered as the definitive temperature 
relative to that derived from the WLP because of the smaller step size and because of a suspected change to the work 
function of the WLP sphere surface. During rapid changes in floating potential, the plasma potential may be outside 
of the sweep range of the NLP, and in this case, no NLP plasma temperature can be obtained. 
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4. Plasma Impedance Probe (PIP) 

The PIP consists of a short dipole antenna electrically isolated from the ISS. The PIP measures the electrical 
impedance (magnitude and phase) of the antenna at 256 frequencies ranging from 100 KHz to 20 MHz in one 
second. Electron density, electron-neutral collision frequency, and magnetic field strength can potentially be 
deduced from these impedance measurements. 

B. Enhanced Capabilities 

In addition to having more instruments and more redundancy than the original FPP, FPMU has greatly improved 
time resolution. Rather than 20 second time resolution on floating potential measurements (in the original FPP 
analysis), FPMU has 128 Hz resolution on floating potential. It was this enhanced time resolution that enabled 
FPMU to detect the rapid charging events that are the focus of this paper. 

C. Discovery of RCEs 

Craven et al. 19 provide a survey of the different types of charging events observed by the FPMU. The normal 
charging events previously observed by the original FPP are also seen by the FPMU at every eclipse exit when the 
PCUs are not in operation. A new type of event, also described by Craven et al. 19 , was dubbed the Rapid Charging 
Event (RCE). The distinguishing characteristic of RCEs, compared to more normal charging events, is the very 
rapid (< 10 seconds) rise time for the charging, which could not have been detected by the original FPP. Typically, 
RCEs only last for less than 30 seconds, and have a faster rise time than decay time. A typical RCE is shown in 
Figure 1 below. The charging is negative, but in the rest of the paper we will refer to its amplitude as a positive 
number. 


2008/235/21:26:17 



In addition to the rapid onset of charging, it is mentioned in Craven et al. 19 that RCEs have a knee in their voltage 
decay profile several seconds after their onset. Examination of the ISS power system telemetry has provided the 
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explanation for this knee. It is ISS power system power management, which starts shunting strings automatically 
within 10 seconds after the onset of current from the arrays at eclipse exit. It is expected that as power demands on 
the arrays go up with time, this shunting will be less common, and the knee may become less prominent. 

In addition, Craven et al. 19 mention the fact that RCEs only appear to occur for very low values of plasma density, 
and with one possible exception, have only been seen when ISS was coming out of eclipse. The one exception, an 
apparent incident at eclipse entrance, has since been discounted, as it is slower than typical RCEs and does not 
adhere to the typical RCE waveform. This event may be part of a separate class of charging events that occurs near 
the terminator. These and other examples of this type of charging event will be addressed in a later paper. 


IV. Rudiments of LEO Charging 

A. Current balance 

Because of the incredibly short electron response time (microseconds at most) a current balance is always 
maintained between ions being collected by some part of the ISS and electrons collected somewhere else. At 
equilibrium, the surfaces of dielectrics are at the potential where thermal electron collection is exactly balanced by 
ram ion collection, which is at most a few kT e away from the plasma potential. For the ionospheric plasma, with 
electron temperature T e of perhaps 0.3 eV at most, this means that insulating surfaces should come to equilibrium at 
most about a volt negative of the plasma. 

Of course, conductors will try to maintain a constant potential across their length, which may or may not be near the 
plasma potential, depending on the impressed potentials of the power system. In particular, when one has a high 
voltage array, such as the 160 V case for ISS, conductors and solar cells on one end of the sunlit array may be at the 
power system voltage different in potential from those on the other end. For ISS and most other satellites, the 
negative end of the power system is grounded to the structure, so the ISS structure may be 1 60 V negative of the 
most positive end of the ISS array. We call the difference of the structure potential from the plasma potential the 
“floating” potential of the structure. It will depend on the relative currents of electrons to the array and ions to the 
array and structure. Most space solar arrays tend to float about 90% negative, because the positive (electron- 
collecting) end of the array collects the lighter electrons about ten times as fast per unit area as the negative (ion- 
collecting) end of the array can collect the heavier ions. 

1. Ram Ion collection 

ISS in its orbit travels faster than the thermal velocity of the (mostly oxygen) ions it flies through. However, it is 
much slower than the thermal velocities of the electrons (the mesosonic condition). In general, this means that only 
the ram ions need be considered in current balance. 

a. Focusing factors (with array and structure geometrical dependences) 

A large flat conducting plate traveling through the ionosphere along its normal will collect the ram ion current 
density (current per unit area). 


J i = N i vq, (2) 

where Nj is the ion number density (usually equal to the electron number density N e ), v is the plate velocity, and q is 
the ion charge (usually one electron charge, e). This can also be expressed as Ji = N e v e. The plate rams into the 
ions with a speed equivalent to 4.89 eV in energy. Because the orbital speed is much greater than the ion thermal 
speed, ions effectively move in straight lines with respect to the velocity vector of the plate. 

A collecting conductor of radius small compared to a Debye length (see the next section below) at a negative 
potential V with respect to the plasma will see focusing of the ions. We define the current density at the surface of 
the collecting conductor as a primed variable. Ion focusing gives us the equation (next page): 
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Ji’ = N e v e (l+|V|/Vo) 1/2 , if the ions are collected by a wire perpendicular to the velocity vector, or 


( 3 ) 


Ji’ =N.ve (1+|V|/V 0 ), if the ions are collected by a sphere. (4) 

See for example, Chen 17 . This can make ion collection much more rapid than without focusing if the ramming 
conductor is surrounded by space or by an insulator with surfaces at nearly the plasma potential. Surfaces in the 
wake will see very little ion collection because of the mesosonic condition. 

2. Thermal Electron collection 

Because the electrons are traveling much faster than the ISS, they can be collected from all angles. At first, this 
makes it seem like they can be effectively collected even in the wake, where the ions are excluded. In actuality, 
however, electrons are also excluded from the wake by the negative space charge that builds up there, so effectively, 
electron collection can only take place on ram or side surfaces. Where the electrons can be collected, their one-sided 
thermal current density per unit area is given by the equation 32 : 

J e = N e e (kTe/27rm e ) l/2 = 2.49 x 1 O' 16 N e (T e ) 1/2 , with T e in Kelvins. (5) 

This is the current density that would be seen by a surface at the plasma potential. However, in actuality, three 
dimensional bodies in plasma are surrounded by a plasma sheath, such that if a thermal electron comes within the 
sheath, it must hit the surface and be collected. For surfaces near the plasma potential, the sheath envelops the 
surface to a distance on the order of a Debye length, given by the following equation 31 : 

X d = (eo kT e /N e e 2 ) 1/2 = 7430 [Te/N e ] 1/2 , with T e in eV. (6) 

If the surface is at some positive potential V, the sheath will be bigger, and in the (orbit-limited) case where the 
sheath radius is much bigger than the radius of the collecting body, the sheath radius is given approximately by the 

. it) 

equation : 


r^ = 0.78 (V/T e ) 1/2 Ad = 5.80 x 10 3 (V/N e ) 1/2 . (7) 

a. Sheath effects (with geometry dependence) 

Hereafter, we use Kelvins for T e . Because all thermal electrons entering the sheath are collected, the effective 


electron thermal current density at the surface of the collecting body will be: 

J e ’ = (Agh/Abody) J e , which for a cylindrical electron collector is (8) 

J e ’ = (r S h/rbody)J e = (2.49 x10' 16 )(5.80x 10 3 )(V N e T e ) 1/2 /r body , or (9) 

J e ’ = 1 .45 x 10' 12 (V N e T e ) 1/2 / rbody , (10) 

and for a spherical electron collector is: 

J e ’ = (WW 2 J e = (2.49x10'’ 6 )(3.3 6x1 0 7 ) V (T e ) 1 /2 /r body 2 (1 1) 

3. Charging times for dielectrics 


Most ISS conductors are covered with a dielectric material, usually for thermal reasons. The conductor underneath 
forms one surface of a (parallel plate) capacitor, with the other surface being the exposed surface of the dielectric. 
Like any capacitor, it takes a certain amount of charge to build up on this surface capacitor before the underlying 
electric field will be neutralized. This charge is Q = C V, where V is the potential of the underlying conductor, and 
C is the capacitance of the material thickness. Dividing by the area, we have Q/A = V C/A, and when charge is 
collected from the plasma, dQ/A dt = J = dV/dt (C/A). Thus, charge builds up so that dV/dt = J/(C/A). For 
dielectric materials of uniform thickness, C/A = e 0 tc/d, where d is the thickness of the dielectric, e 0 is the permittivity 
of free space, and k is the dielectric constant of the material. Taking values which we will later see are relatively 
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typical for the RCEs (k = 5, a -22 V charging voltage, a plasma density of 2x1 0 10 and an electron temperature of 
2000 K), we have for the various ISS charging surfaces the values in Table 1 below: 

Table 1. Timescales for ISS charging in a plasma of N e = 2x1 0 10 and T e = 2000 


Type of Surface 

d (dielectric 
thickness, m) 


Type of charging 

1/t (1/s) 

t (sec) 







ISS Structure type 1 

1.30E-06 

3.40E-05 

Ram Ion 

3.29E-02 

3.04E+01 

ISS Structure type 2 

1.30E-05 

3.40E-06 

Ram Ion 

3.29E-01 

3.04E+00 

ISS Solar Array 

1.50E-04 

2.95E-07 

Ram Ion 

3.80E+00 

2.63E-01 

ISS Structure type 1 

1.30E-06 

3.40E-05 

Focused Ram Ion 

2.30E-01 

4.34E+00 

ISS Structure type 2 

1.30E-05 

3.40E-06 

Focused Ram Ion 

2.30E+00 

4.34E-01 

ISS beta cloth 

7.60E-04 

5.82E-08 

Ram Ion 

1.92E+01 

5.20E-02 







ISS CIC (cell gap) 

7.60E-05 

5.82E-07 

Thermal Electron 

4.32E+00 

2.31E-01 

ISS Solar Array 

1.50E-04 

2.95E-07 

Thermal Electron 

8.53E+00 

1.17E-01 


4. Magnetic field effects 

ISS in its orbit passes through the Earth’s magnetic field. The magnetic field can constrain the lightweight electrons 
more easily than the heavier ions. In our initial attempt at understanding the RCEs, we ignore the effects of the 
magnetic field on current collection. 


V. Theory for RCEs 

A. Conditions for Success 

The attributes of RCEs stated in Craven et al. 19 set the criteria for a successful theory; namely, the amount of peak 
charging, the steep inverse dependence of the peak value on the electron density, the rapid timescale(s) for charging 
and charge decay, and the occurrence of RCEs only at eclipse exit. 

1. Early “sheath” charging (ps timescales) 

When ISS is in eclipse, and the solar arrays are not producing a distributed voltage, all surfaces of ISS are 
effectively negative by a few times the plasma potential, at most a volt or so. As ISS comes out of eclipse, at some 
insolation value the solar array voltages can switch on very suddenly, placing the positive and negative extrema of 
ISS about 1 60 volts apart. In what follows immediately, we will describe what would happen if the array turn-on 
time were instantaneous. 

First of all, because the dielectric surfaces are capacitively coupled to the underlying conductors on short timescales, 
each of the solar cell coverglass surfaces will take on the same potential as the underlying cell. This means that 
before dielectric surface charging can take place, all surfaces will be able to collect ions or electrons as if they were 
conductors. In this condition, the ISS structure would float highly negative, because the solar arrays, with their 
much larger surface areas, would dominate the charging. And remember, if both the ion and electron collection of 
an array are unhindered, it will float about 90% negative. This would put the ISS structure at about 140 V negative 
at eclipse exit. 

Remember, too, that the capacitance of ISS to space is very small (effectively being the capacitance of a vacuum 
layer the thickness of the sheath around ISS). Both ion and electron collection can charge up the sheath capacitor in 
less than 10 microseconds, so for an instantaneous array turn-on, ISS would charge to about -140 V within 10 
microseconds. This has not been observed. Observations show that the voltage turn-on time for the RCEs is about 
0.1 second. Why is this? 

It has been found that for an individual solar cell, the voltage turn-on time is very short (probably less than a 
millisecond). For a distributed array of cells, however, differences in individual cell activation insolation values 
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may spread out the full array voltage turn-on time. From simple geometry of the ISS orbit and the angular diameter 
of the sun, one can find that the time for the array to reach lull insolation is at least 8.0 seconds and can be 10 times 
that long. Thus, if there is a distribution in illumination for activation of individual cells of one part in 80, it would 
have the effect of smearing out the other effects we mentioned above over 0.1 second, by which time other charging 
effects will have become important. 

In what follows, we assume that the array voltage turn-on time is greater than about 1/10 second (which agrees with 
the observed sharpness of the voltage turn-on time), so the initial very rapid structure charging (sometimes called 
absolute charging) does not develop. We note, also, that ram ion charging of the ISS beta cloth happens within 
about the first 1/10 second, so that for conductors surrounded by betacloth, ion focusing will occur from the 
beginning. 

2. Special initial condition (uncharged coverglasses) 

Thus, we will take as zero the time of fully developed voltage on the array, which we take to coincide with the time 
it takes for ion focusing to be fully developed on ISS betacloth-surrounded structure conductors. The surfaces that 
have not been charged to near plasma potential are the coverglass and Kapton®-covered parts of the solar arrays and 
all structures coated with type one and two dielectrics (anodized aluminum). 

a. Initial current balance equation 

Ion collection to ground will happen on exposed conductors, of which there are three main types: 

i. Structural conductors surrounded by betacloth. Although the total area of conductors on ISS is 
estimated at 20-35 m 2 , it is not known what fraction is composed of small conductors where focusing 
will be important. Because small conducting areas are surrounded by betacloth at about the plasma 
potential, we assume spherical ram ion collection on the forward hemisphere only, and Jj’ = (l+|V|/4.89) 
N e v e/ 2 = 1.2 x 10‘ 15 (l+|V|/4.89) N/2 . At a density N e of 2 x 10 10 , Jj’ = 2.4 x 10 s (l+|V|/4.89)/2 
amp/m 2 , and if we have charging of -22 volts, say, this amounts to 8x1 O' 5 amp/m 2 . 

Ion collection will be enhanced by photoemission on the conductor surfaces. This amounts to about 2 x 
10" 5 amp/m 2 in full solar ultraviolet. Not a trivial amount, it should be added to the ion collection of 
structure elements when in full sunlight. However, ultraviolet sunrise on ISS occurs many seconds later 
than infrared sunrise, which is responsible for turning on the solar array voltages. (The spectral 
response 18 of ISS solar cells peaks in the wavelength range of 900 to 1000 nm, the near-infrared.) This 
delay can be seen in Figure 2, below. 



Figure 2. Geometry for atmospheric absorption. 
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( 12 ) 


COS#! 


cos # 2 


r E + ALT 


r*+d 


Atm 


r v + ALT 


(13) 


Here, we have assumed an orbital inclination of 0 degrees. For greater inclinations, the effects will be 
greater. The vertical line is perpendicular to the solar direction. ALT is the ISS orbital altitude. r E is 
the Earth’s radius. ISS first sees the sun through the entire clear atmosphere (and at infrared 
wavelengths) at 0\. For an upper limit of complete ultraviolet absorption of d atm , ISS first sees the 
ultraviolet sun at # 2 . The orbital time between d\ and # 2 is the time delay for the ultraviolet 
photoemission after the first voltage turn-on. 


Putting in ALT = 350 km, r E = 6378 km, and a conservative estimate for d atm of 30 km (the lowest level 
of the ozone layer), we find that # 2 - d\ = 0.85°, and the time delay is about 13 seconds of time. This is 
beyond the peak of all RCEs treated in this paper, and so we can safely assume that photoemission by 
ultraviolet light is unimportant for RCEs at a clear atmosphere eclipse exit, because the ultraviolet starts 
reaching ISS long after the infrared emission has initiated the RCE. If the atmosphere is cloudy at the 
eclipse exit, the photoemission and voltage terminator may come at the same time, the relative 
contribution of photoemission may be greater by the amount indicated, and this would make for a lower 
RCE voltage peak. 


For spherical collection at a cloudy eclipse exit, the local effective ion collection current density on a 
small conductor surrounded by an insulator is: 

Ji’ = [1.2 x 10 15 (l+|V|/4.89) N/2 + 2 x 10‘ 5 ] A/m 2 . (14) 

For spherical collection at a clear eclipse exit, 

Ji’ = [1 .2 x 10 15 (l+|V|/4.89) N e /2] A/m 2 . (15) 

In the following, we will assume clear eclipse exits, since we are interested in the maximum amount of 
charging that may occur (and thus the smallest ion current density). 

ii. Wires on the solar array masts. Each solar array mast has over 500 m of tensioning wires or rods 42 . 
Assuming a 2.38 mm cross-section for the wires, this gives a total collecting area of about 7.5 square 
meters for the 6 masts with three arrays deployed. However, geometrical factors will limit their average 
collection to about 14 that of wires all oriented perpendicular to the ram. Thus, these will collect ions 
like cylinders, according to the equation: 

Ji” = N e v e (l+|V|/V 0 ) 1/2 /2 = 1.2xl0' 15 (l+V/4.89) I/z N e /2 (3a) 

iii. The edges of the solar cells on the arrays. Before the ram ions have had a chance to neutralize the 
array coverglasses, this will be unfocused ram ion collection, with J ; = 1.2 x 10' 15 N e . For sake of 
comparison with other current densities, we can put in N e = 2 x 10 10 to find an unfocused ram ion 
current density of about 2.4 x 10' 5 amps per m 2 , or 2.4 nA/cm 2 . 

As we will see later, the electron current density on the cell edges far overwhelms that of the ions, even 
for focused ions, and so ion collection on the array will be unimportant. For completeness, however, we 
include unfocused ram ion collection on the array in the initial charging equation. 

Electron collection will only occur on the solar cell edges, and will, in general, be collected from a 
sheath. If the coverglasses are completely un-neutralized (as at t = 0), the solar cells can collect current 
at all of their edges. Each string on the array has 400 cells, each of which has about 0.32 m of cell edge 
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of 6 mils uncovered thickness = 1.5 x 10' 4 m, for an area per cell of 4.8xl0' 5 m 2 and a string edge area of 
1 .9 x 10‘ 2 m. There are 82 strings per wing, and we assume four wings or two arrays, so the total edge 
area is 6.3 m 2 . When the coverglasses are completely neutralized, the electron collection will be almost 
completely choked off, and only the outer edges of the outermost cells will be electron collectors. This 
restricts the total effective edge area by a factor of 8. Thus, the geometry of the ISS solar array is 
extremely important, and is taken from Kerslake and Scheimann 18 and presented in Figure 3 below. 


Factory Verification - Wing Assembly (Typical Circuit) 
400 Series Connected Solar Cells 



Figure 3. Geometry of the ISS Solar Array Wing Assembly. 

As can be seen here, each 40 cell module (measuring about 32 cm x 80 cm), is connected in series with 
its neighbors in an across-the-u, back-the-u fashion, with the legs of the u separated by almost one cell 
width. For a structure charging value of 35 V, one end of the sunlit array will be about 125 V positive, 
leading to a sheath radius of 


r sh = 5.80 x 10 3 (V/N e ) 1/2 = 0.46 meters at N e = 2xl0 10 . (16) 

Thus, we expect the sheath, at its broadest point, to be wider than one leg (or both legs) of the u, and so 
the local sheath will be a tube shape surrounding the center of one leg of the u. At V = 125-80 = 45 V, 
the turn back of the u will have a sheath radius of about 0.28 m, so the sheath diameter is still larger 
than the width of the u leg. We feel justified in taking the electron current collection to be from a sheath 
that is a circular cylinder of varying radius along the array length. So, locally we have thermal current 
collection on the surface of a half-cylinder of radius r sh , which then is focused onto a flat plate of width 
0.32 m. J e ’ is then the ratio of these two areas, or 


So, J e ’ = Je r sh (7t/0.32), or 

(17) 

J e ’= 1.42 x 1 O' 11 (V N e T e ) 1/2 . 

(18) 


Whatever the value of the negative structure charging, which we give the value cp, the fraction of the 
array that is positively charged (and electron collecting) will be (1 60-cp)/l 60 and the fraction that is 


negative (and ion collecting) will be cp/160. 

We are now ready to write the current balance equation. Initially, 

JiA SA (cp/160) + Jj’A ISS + Ji”A mast = J e ’A SA (160-cp)/160 = J/A sa - J e ’A SA (cp/l 60), (19) 

where A mast is the mast wire area. 

A sa (cp/1 60)(Ji + J e ’) = J e ’A SA - Jj’Ajss - Ji”A mast , so that (20) 

Cp = 160 [(J e ’A SA - Ji’Arss - Ji”A mait )/(Jj + J e ’)]/A SA , and (21) 

cp = 160 [J e 7(Ji + J e ’)] [A sa - (J,7J e ’)A iSS - (J,”/J e ’)A mast ]/A SA . (22) 
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Here, J e ’ is the average value over the electron collecting part of an array. For all realistic cases, J e ’ » 
Jj,, so the charging equation becomes 

q> = 160 [1 - (Jj’Ajss + Ji”A mait )/ J e ’A SA ]. (23) 

As a sanity check, this guarantees that for a total ion current greater than the electron current, there will 
be no charging. 

We can calculate the average J e ’ on the solar array as 


J e ’ = 1 .42 X 10- 1 1 (Ne T e ) 1/2 (V ,/2 avg ). (24) 

Doing a numerical integration of W v \ we find that (V y2 av g)~ (0.943)(V avg ) 1/2 , so that 

J e ’ = 1.34 x 10 " (N e T e ) 1/2 [( 1 60-cp)/2] 1/2 = 9.48 x 10‘ 12 (N e T e ) 1/2 (160-cp) 1/2 (25) 

And finally, from (23), 

(p/160 = [1 - (Ji’A IS s + Ji”A mast )/ Je’A SA ], so (23) 

(p/160 = 1 - [(N/T e ) 1/2 / (9.48x1 0 12 )( 1 60-(p) 1/2 ] [(J,’)A ISS + (J,”)A mast ]/A SA , and (26) 


(p/160 = 1 - [(Ng/T e ) / (160-(p)] l/2 [1.30 x 10 4 (1 + (p/4.89) (A ISS /2A sa ) + 1.30 x 10 4 (1 + cp/4.89) 1/2 (A mas /2A SA )], 

(27) 

or, 

(160-(p)/160 = [(N/T c )/ (1 60-<p)] l/2 [1.30 x 10' 4 ] [(1 + (p/4.89) (A^A*) + (1 + cp/4.89) 1/2 (A mast /2A SA )]. (28) 
Finally, 

( 1 60-(p) 3/2 = 160 (Ng/Te) 1/2 [1.30 X 1 0 4 ] [(1 + (p/4.89) (A ISS /2A sa ) + (1 + q>/4.89) 1/2 (A ma5 /2A SA )]. (29) 

Here, we can solve for (p = 0 and see that there will be no charging for 

(Ng/T e ) 1/2 > 9.73xl0 4 [2 A sa /(A iss + A mast )]. (30) 

From eqn. 29, we can see that (p depends only on the areas and (Ne/T e ) 1/2 , a fact that will come in handy 
later when we compare theory with observation. 

Going forward, we will assume that for the time period we have chosen (see VI. below), with two fully 
functioning arrays and one not tracking, A IS s = 20 m 2 , A mast = 7.5 m 2 , and A sa = 6.3 m 2 . With these 
assumptions, an excellent numerical approximation to the complicated charging equation (29) is the 
following : 


In (p = 5.194 - 0.02415 V[(N/r e ) 1/2 ] 


(31) 


Putting in N e = 2 x 10 10 and T e = 2000, typical parameters for RCEs, we solve for (p = 46 V. This is 
similar to the charging values seen for some RCEs. We emphasize here that this number is extremely 
sensitive to the total conducting area of ISS as well as that of the solar arrays, and is thus uncertain by the 
uncertainty in those quantities. 

It is important to note that the charging equation predicts a steep dependence on (N/r e ). For instance, 
putting in a high value of N e (3 x 10 10 ) and a low value of T e (1000) compared to typical RCEs, we find (p 
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= 30 volts. Conversely, for N e = 1 x 10 10 and T e = 3000, we find cp = 64 volts. This is the behavior first 
noted 19 for RCEs in peak voltage versus N e . 

As the coverglasses charge, the electrons will be focused down onto the cell gaps, and they will also 
charge, effectively cutting off electron collection to interior cell edges as noted above. In addition, as ions 
that are being focused on ISS surfaces but miss the conductors overcharge the adjacent areas, this will 
choke off ion collection by the structure. Thus, both electron and ion collection will decay with time. 
From equation (21) it is clear that it is the difference between the electron collection and the ion 
collection that determines the charging. At t = 0 there is no charging, so they must be equal. After a long 
time has elapsed, again the charging is small, so both are equal. Since both electron and ion collection 
decay with time, we expect the charging curve to be a difference of two exponential decays of the same 
amplitude but different timescales. The (p found above will then be the amplitude from which the two 
decays start out. That is, 

V obs = -cp [e' at - e' bt ], (32) 

Here, we expect the parameter “a” to be related to the decay of ion collection, since according to Table 1, 
ion-collecting dielectrics on surface structures will take a longer time to charge up than will the solar 
array coverglasses. Thus, we expect “b” to be related to the decay of electron collection due to coverglass 
charging. 

b. Dependencies of charging amount 

• Density and Temperature dependence. According to the charging equation, the amplitude of 
charging should depend almost entirely on the collecting areas and (N/T e ) 1/2 . 

• Time dependence. The charging peak should be the result of the difference of two exponentials as 
discussed above. 

• Amplitude versus observed amount. 

Vobs = - «p [e' at - e‘ bt ]. We define the ratio a/b = x. It is cp and x that determine the shape of the charging 
peak, aside from the time as a scaling factor. 

• Dependence on array angle. 

Since the beta angle of the ISS orbit varies from ~ -70° to +70°, a tracking array will not always be 
pointing directly into the ram at eclipse exit. The amount of electron collection depends on the cosine 
of angle the normal of the array makes with the velocity vector, so there should be a cosine dependence 
of cp on array ram angle at eclipse exit. 

• Time of peak voltage. 

Differentiating eqn. (32) with respect to time and setting the derivative equal to zero, we find that the 
peak should occur at a time 


tp k = (In a - In b)/(a-b). 


(33) 


• Rate of change of voltage. dV obs /dt = -cp [-ae' at + be' bt ], and differentiating this with respect to time 
and setting it equal to zero, we find that the maximum in d|V ob8 |/dt is at time 

t = 2 (In a - In b)/(a-b) = 2 1*. (34) 

It should be noted that there is also a maximum at a time immediately after t = 0 (to in Figure 1). 
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c. Transition to “normal” charging - final current balance after coverglass charging 

“Normal” charging for ISS has electron collection severely restricted by charging of the array coverglasses, 
especially for electron energies above about 0.2 eV. Thus, it is in general lower than RCE charging and 
goes up exponentially with lower electron temperatures 12 . However, this will be the charging level to 
which the RCEs decay. 


VI. RCE Observations and Comparison to Theory 

The FPMU has been operated intermittently on ISS since mid 2006. Since that time, it has amassed a fairly large 
statistical base of RCEs, along with ancillary data (N e , T e , peak voltage, beta angles, etc.). In what follows, we will 
analyze RCEs from the period Jan. 28 through Nov. 4 of 2008. In that time period, about 80 RCEs were observed. 
Many of the RCEs had complex shapes in time. The difficulty in determining an all encompassing theory for them 
is the many variables involved that must contribute to the variety of complex shapes observed. However, there is a 
certain “canonical” shape that is apparent and in this paper we will concentrate on the events with that particular 
shape (see Figure 1). For our analysis we are using the N e and T e values averaged over a 5 second time period, just 
before eclipse exit. We assume that for the RCEs of simple shape, these parameters do not vary during the time 
extent of the RCE. We have adopted this strategy because (as was mentioned before) our best determinations of Te 
come from the NLP, which is unfortunately not able to keep up with the rapidly changing floating potential during 
an RCE. Thus, we wish to determine that a simple model can be constructed for RCEs of simple shape (where 
conditions are not rapidly changing). In the final analysis, a uniform sample of 42 RCEs of the simple “canonical” 
shape were chosen. Thus, our sample size is 42 events. In the following, the units of T e are Kelvins. 


A. RCE Shapes - fits by two exponentials 

A curve fitting routine 20 (TableCurve 2D v4.07) was used to fit various peaked functions to the RCE waveforms, 
using data with a time resolution of 128 values per second of time. Since the waveforms decayed down to a non- 
zero baseline, it was decided to fit only the first 10 to 20 seconds of the waveform, to prevent contamination from 
“normal” charging. In all cases, almost all of the rising voltage and a good deal of the falling voltage regions were 
fit. Of the 37 peaked functions available in TableCurve 2D, the so-called Intermediate Peak Function (a difference 
of two exponentials) gave the best fit to almost all of the RCE shapes for the 5 free parameters available (amplitude, 
two exponential time constants, zero offset, and zero time). All RCEs analyzed had correlation coefficients in the 
range from 0.895 < r 2 < 0.998. Two of the representative data fits are shown in Figures 4a and 4b. In both these 
plots, time in seconds is on the abscissa and negative voltage in volts is on the ordinate. 



Figure 4a. A representative RCE and TableCurve fit by the difference between two exponentials. 


14 

American Institute of Aeronautics and Astronautics 





Figure 4b. A representative RCE and TableCurve fit by the difference between two exponentials. 

From these fits, values for (p, voltage offset, a, and b were obtained for our sample of 42 RCEs. 

1. Values and ratios. 

The voltage offset found from the model fits was highly correlated with the approximate estimate found by eye. 

Values of the parameter “a” from our fits ranged from 0.43 to 5.5 s' 1 , (for 1/e times of 0.18 to 2.3 seconds) with a 
mean of 1.61 and median of 1.24 s" 1 (0.62 and 0.80 seconds, respectively). It is instructive to inquire as to the N c 
and T e dependences of “a” from our fits, so as to shed light on the physical mechanism for this decay. For each 
RCE, the density and temperature were determined as the average of 5 FPMU measurements from just before 
eclipse exit. In Figure 5 below is the relationship of “a” to N e for our sample. For a sample size of 40, the 0.0005 
confidence level is at r = 0.50 (r 2 = 0.25), so for our sample of 42 we are at a confidence level 41 of better than 0.0005 
for “a” vs. N e . 

a(s _1 )vs Ne 



NeflO 11 m- 3 ) 

Figure 5. The value of “a” versus electron density. 
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The dependence of “a” on T e , on the other hand, is not statistically significant at even the .05 level. Since “a” is 
strongly correlated with N e , “a” may be the decay time constant for ion collection, as we reasoned before. 

Values of “b” from our model fits ranged from 0.056 to 0.60 s' 1 (for 1/e times of 1.7 to 18 seconds), with a mean of 
0.21 and median of 0.18 s' 1 (4.8 and 5.6 seconds, respectively), “b” shows significant correlations with both N e and 
T e , and for simple combinations ofN e and T e , the best correlation found was an exponential of “b” with (N e T e ) 1/2 . 
See Fig. 6. 


b (s-l)vs V (Ne*Te) 



V (Ne*Te) 

Figure 6. The value of b versus V(NeTe). 

Since electron collection is associated with (N e T e ) 1/2 , “b” may thus be the time constant for decay of electron 
collection, as we had also reasoned before. 

Values of the ratio of a/b ranged from 2.40 to 52.7, with a mean of 10.6 and median of 6.3. 

We will call the value of (p obtained from the TableCurve fit to each RCE “cpm”- <pnt ranged from 7.35 to 91 .9 volts, 
with a mean of 28.7 and median of 23.2. 

2. Density and temperature dependences 

The amplitude <pf lt from the TableCurve fits is related to (N</T e ) 1/2 as shown in Figure 7. 
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q> vs V (Ne/Te) 



— Phi 
a Phi(f it) 

Linear {Phi(fit)) 


= 0.3309 


V (Ne/Te) 

Figure 7. Amplitude versus V(Ne/Te) 

In this figure Phi is the value of amplitude cp predicted by Eqn. (31) which closely approximates the charging 
condition given by Eqn. (29). Here it can be seen that the theory predicts the peak value of the amplitude versus 
(N/Ts.) 172 very well (except for two discrepant points), and also that the amplitude of the RCEs is linearly correlated 
with (Ne/T e ) l/2 with a highly significant correlation coefficient of r 2 = 0.3309. 

We will discuss below why not all amplitudes of the RCEs lie along the predicted line. 


3. Dependence on ram angle of arrays 

Because the solar arrays are the sole electron collectors, and they are not oriented directly into the ram direction, 
their efficiency in producing charging should depend on the ram angle of the solar array at eclipse exit. In Figure 9, 
we see the ram angle dependence of the amplitude. Here, we have normalized the observed amplitude by dividing 
by <p from equation (31) above, so as to take out all other dependences. 
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<pf it /<p vs Solar Array Ram Angle 



♦ Phi (fityPhi 
— 1.35 cos(ram) 


Also plotted here is a cosine function, normalized to fit the maximum values of amplitude. It can be seen that the 
value of the peak amplitudes does indeed fall off as the cosine of the ram angle, as one would predict from the 
geometry. As in figure 7, the same two discrepant points spoil the excellent correlation. 

4. Other factors - ISS orientation, etc. 

There are many possible reasons why not all of the RCEs take on the amplitude calculated from the theory. First of 
all, the eclipse exit for a given RCE may have been under cloudy conditions, where the onset of voltage was 
coincident with the onset of ultraviolet photoemission. In this case, as we have noted before, the absolute amount of 
charging would be lowered by the negative charge loss due to photoemission. 

Secondly, there may have been changes of T e and N e during the RCE duration. As ISS comes into daylight, one 
naturally expects N e and T e both to increase. Since ISS moves about 1 orbital degree in 1 5 seconds, in order for the 
change to be important during the 20 seconds or so of the RCE, it would have to happen during about the first one or 
two degrees after eclipse exit. This is possible, if the gradient of T e and/or N e are steep enough at the terminator. 
We may be seeing the effects of high gradients in N e and T e in the RCEs with complex structure. Also, so-called 
“density dropouts” have been observed in the ionosphere, where the density decreases and the temperature increases 
(so as to equalize the pressure) on very short distance scales. Examination of the time history of T e and N e during 
some of our observed RCEs shows that this effect may make for an occasional discrepant point. 

Thirdly, during the time interval over which the data for the model fits was obtained, there were really three 
functioning arrays on ISS, but the third array was not tracking (and, in fact one half was pointing in a different 
direction than the other). Thus, depending on beta angle, it could supply more or less electron collection from the 
illuminated half or ion collection area from the unilluminated half. Also, drag reduction modes of flying ISS, 
including the so-called night-glider mode, with the solar array wings partially or fully feathered to the velocity 
vector, can cause large variations in the electron collection at eclipse exit. 
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Fourth, power management may have taken place during the RCE. Usually, at eclipse exit, all strings on the active 
arrays are unshunted (active). However, shortly after the arrays start producing power, power management 
algorithms start turning off strings to maintain the only the current needed to keep the batteries fully charged. 
Depending on the state of battery charge, this may occur sooner or later in the orbit. Shunted strings are 
disconnected from system ground, and will not contribute to ISS current balance. 

Fifth, small changes in ISS orientation at eclipse exit (there is a deadband in angle around the nominal orientation, 
with corrections at the ends of the deadband) may make a big difference in ion collection. Turning the ion collection 
areas of the Russian segments even a little bit into or out of the ram can be important. Also, during reboost and at 
other times, the ISS orientation can vary from the normal +XW attitude. 

Sixth, in our simple-minded derivation of the ion and electron collection, we have ignored all magnetic field effects, 
including vxEM effects. Depending on the magnetic latitude and longitude of eclipse exit, this may change the 
relative proportions of array electron and ion collection. 

Some of these effects may be estimated from detailed examination of the ODRC data, but that must remain for 
forward work. Suffice it to say that the simple theory presented here works well for the RCEs of simple shape, and 
there is enough flexibility in orbital conditions at eclipse exit to make it plausible for all of the RCEs. 


VII. Predictions for ISS Final Assembly and Peak Solar Activity 

As the number of active and tracking electron-collecting arrays on ISS increases, without a concomitant increase in 
the structure ion collecting area, we may expect the RCEs to be of higher amplitude. Putting in values for the 
increased electron collection and ion collection area (A sa = 12.6 m 2 , A mast = 10 m 2 ) into the calculation of cp above, 
and using a numerical approximation, we have the equation below and Figure 9 showing predictions for RCE 
amplitude versus (N,/r e ) 1/2 for the 4 array case: 

In cp’ = 5.2147 - 0.01741 V[(Ne/T e ) 1/2 ] (35) 

(p vs V (Ne/Te) 



Figure 9. Peak amplitude predictions for 4 active arrays. 


♦ Phi (fit) 

— Phi' 4 arrays 
— Phi 2 arrays 
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Figure 9 shows that for the full 4 arrays in operation, peak RCEs of amplitude at least 81 volts are predicted, with 
routine values in the 60-70 volt range. Putting in average values of “a” and “b” into equation (4), and calculating a 
peak voltage at that time, these amplitudes predict that the peak voltages are likely to be in the 45-55 volt range. 
Note also that for 4 solar arrays active, eqn. (30) shows that RCEs extend to higher values of (Ne/T e ) than for 2 
arrays, so we may expect that in 4 solar array operational conditions, there will be a higher percentage of eclipse 
exits with some degree of rapid charging. And, from Figure 8, we might expect that for a maximum ram condition 
the peaks will be about a factor of 1 .3 times as great as from cp alone, and this makes for a maximum peak amplitude 
of about 105 volts, with peak voltages in the 75 volt range. 

These predictions are based on low solar activity, when it can be expected that (N</r e ) is low when coming out of 
eclipse. At times of high solar activity (Ne/T e ) will be higher, on the average, and it is less likely that severe RCEs 
will occur. However, as has already been pointed out 21 , during solar geomagnetic substorms (which occur 
preferentially during times of high solar activity), the variability in both N e and T e are greatly enhanced, and periods 
of low (Ne/T e ) may still be expected. 

A. Probability of dielectric breakdown from RCEs 

Because RCEs represent events where the structure potential becomes highly negative before all dielectric surfaces 
are neutralized, we should not expect that the electric fields developed across dielectrics are as great (compared to 
the amplitude of structure charging) as they would be for normal charging. It is the potential difference between the 
underlying structure and the dielectric surface that produces the electric fields that are important in dielectric 
breakdown, after all. One may be justified in asking, then, what is the potential difference for dielectric materials 
dining an RCE? In general, dielectric surfaces will be neutralized by ion collection, so one may compare the 
instantaneous structure potential with the surface potential on the dielectric, which is changing from cp to near zero 
on the dielectric ram ion charging timescale for that particular dielectric. Looking at the time equation for charging, 
we can postulate that 


V P k - Vdi = <P [e' at - e" bt - e~ ct ], (36) 

where “c” is the charging timescale for the dielectric (di) in question, and t is the time of the peak from equation 
(33), 


tp k = (In a - In b)/(a-b). (33) 

Identifying the timescale a as the ram ion charging timescale for at least one type of ISS structure dielectric, we find 
for that dielectric, c = a and, 


|V pk - V<ji| = cp [e' bt ], so that at tp k , we have (37) 

|V pk -V di | = cp[(b/a) (b/[a - bl) ]. (38) 

Putting in values typical of a andb for RCEs (1.61 and 0.21, respectively), we find 

|V pk - Vdi| = cp [ (0.13) 015 ] = cp [0.74], (39) 

Thus, we may expect this dielectric to have a peak potential difference across it of about 0.74 of the maximum peak 
amplitude cp. For a 105 volt amplitude, which may sometimes be the case for 4 solar arrays operational, this implies 
a dielectric potential difference of perhaps 78 volts. For ISS chromic acid anodized aluminum (our type 1) material, 
a dielectric breakdown strength of 66-82 V has been found 22,23 . Thus, it seems theoretically possible that ISS type 1 
anodization may break down due to RCEs when there are 4 solar arrays active. However, Schneider et al. also used 
a 10 minute waiting time to detect arcs near the threshold, whereas our charging peaks are on the order of 4 seconds 
wide. If the arcing is purely stochastic, within a 4 second time period we might expect a probability of one in ~1 50 
that any given RCE would produce an arc. If an arc does occur, however, it may be of high amplitude (perhaps 10 
amps) and duration (5 ms or more) and may thus expend several Joules of energy 23 . Type 2 anodization, much 
thicker and more robust, should arc very seldom if at all. We predict, therefore, that after several hundred RCEs on 
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an ISS with 4 active solar arrays, there may occur one or more arcs on ISS type 1 anodization. If the ISS program is 
concerned about the effects of individual dielectric breakdowns, a mitigation technique (either the use of PCUs or an 
alternative such as feathering the arrays 24 ) can be adopted for the first few minutes of eclipse exit, to prevent 
possible arcing from large RCEs or normal charging with 4 active arrays. If arcs do occur, they will be easily seen 
in the floating potentials of a continuously operating high time resolution instrument such as the FPMU. 


VIII. EWB Model 

An EWB model was produced of ISS at Assembly Complete for eclipse exit with the following parameters: N e = 
5.1xl0 10 m‘ 3 , T e = 0.167 eV, ram angle at eclipse exit = 26.3°, latitude of eclipse exit = - 48.5°, altitude = 341 km, 
and v = 7.7 km/s. The solar array electron and ion current collection curve in EWB was adjusted so that it matched 
that expected from the RCE Theory at the plasma densities and temperatures given. Ram ion current collection 
areas were also adjusted to give the correct ratio of ion collection currents to array electron collection currents. The 
EWB model found extrema in the ISS potentials of -1 1 and -52 volts, with the average structure potential of -31 V, 
approximately the value predicted from the RCE Theory. A plot of the potentials is shown in Figure 1 0. The major 
differences in potential across the ISS structure are due to vxB-1. 



Figure 10. Structure potentials on ISS at Assembly Complete configuration for representative plasma 

conditions at eclipse exit, showing vxB*l effects. 
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IX. Conclusions and Predictions 


A. No Ad-hoc Physics 

We have shown that it is possible to model the RCEs with a simple physical model assuming the initial condition of 
uncharged coverglasses, focused ram ion collection and array electron collection through a sheath. 

B. Model Agreement with FPMU Observations 

Our model agrees with the observed maximum amplitude and time behavior of RCEs, as well as their dependences 
on electron density N e and temperature T e . The model also predicts that RCEs will decay into normal charging, 
which is also in agreement with observations. Finally, the model shows that RCEs of large amplitude will not occur 
for large values of (Ne/T e ), also in accord with observations. 

C. Predictions for Dielectric Breakdown at Assembly Complete 

We have calculated whether dielectric breakdown can occur on ISS anodized aluminum surfaces, and have found 
that with 4 active arrays, dielectric breakdown is likely to occur after several hundred RCEs have happened. In 
addition, we find that the probability of RCE occurrence is enhanced by having 4 solar arrays active. 

D. Predictions for Solar Max 

Because (Ne/T e ) is, on average, higher at eclipse exit near the maximum of the solar activity cycle, this should make 
the RCEs rarer than at the present time. Elowever, during solar geomagnetic substorms, the likelihood of RCEs 
increases, due to the wider than normal deviations from average (N/r e ) conditions at those times 21 . 


E. Future Use of Plasma Experiments on ISS 

ISS has shown that it is a wonderful platform for plasma experiments in LEO. Even with instruments that worked 
for only a short time (FPP), or for limited time periods (FPMU), and with a smaller number of active arrays, the 
knowledge that has been gained about LEO spacecraft charging has been dramatic. With the coming of plasma 
instruments on the Japanese Experiment Module and the increased power requirements on the solar arrays, which 
will prevent their shunting after the first few seconds of power generation, the opportunities for investigating ISS 
charging will greatly increase. In order to fully understand the RCEs, the following data and experiments would be 
desirable: infrared and ultraviolet fluxes to help determine the instants when solar array voltages are tumed-on and 
photoemission starts, occasional periodic on/off PCU cycles just after eclipse exit to determine the RCE and 
normal charging baselines, continuous plasma parameter and floating potential determination so we won’t miss any 
RCEs or discharges due to RCEs, and floating potential determinations during docking and periods when Shuttle or 
CEV are docked to see how the presence of the current collection of these bodies changes ISS charging. Because of 
its size, high voltage arrays, and interesting structure-conductor and solar cell geometries, the ISS and its plasma 
interactions will undoubtedly continue to surprise and enlighten us. 
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Typical Rapid Charging Event 



Rapid Charging Event (RCE) Scenario 

When ISS comes out of eclipse, voltage turns on 
rapidly 

Solar cell coverglasses are initially at potential of 
underlying cell because of capacitive coupling 

Electron current choking on solar arrays is not in 
effect, since the coverglasses and gap kapton are 
not near plasma potential 

Until coverglasses can charge, high values of 
current collection allow ISS to charge rapidly to 
high negative voltages 

otential will decay as electrons and ions can 
charge surfaces to near plasma potential 



Timescales for ISS charging in a plasma of 
Ne = 2xl0 10 nr 3 and Te = 2000 K 


Dielectric 

thickness, 


Type of Surface 

m) 

C/A (F/m 2 ) 

Type of charging 

1/t (1/s) 

t (sec) 

ISS Structure type 1 

1.30E-06 

3.40E-05 

Ram Ion 

3.29E-02 

3.04E+01 

ISS Structure type 2 

1.30E-05 

3.40E-06 

Ram Ion 

3.29E-01 

3.04E+00 

ISS Solar Array 

1.50E-04 

2.95E-07 

Ram Ion 

3.80E+00 

2.63E-01 

ISS Structure type 1 

1.30E-06 

3.40E-05 

Focused Ram Ion 

2.30E-01 

4.34E+00 

ISS Structure type 2 

1.30E-05 

3.40E-06 

Focused Ram Ion 

2.30E+00 

4.34E-01 

ISS beta cloth 

7.60E-04 

5.82E-08 

Ram Ion 

1.92E+01 

5.20E-02 

ISS CIC (cell gap) 

7.60E-05 

5.82E-07 

Thermal Electron 

4.32E+00 

2.31E-01 

ISS Solar Array 

1.50E-04 

2.95E-07 

Thermal Electron 

8.53E+00 

1.17E-01 


Geometry of the ISS solar array 
Wing Assembly 




Current Balance Theory - <|> = Amplitude 
(no coverglass charging) 

cp /160 = [1 - (Jj’Ajss + J”A mast )/ J e ’A SA ], 

where 

Je’ = 9-03 X io- 12 (V avg N e T e ) l/2 

(array electron thermal collection from sheath) 

]■ = [1.2 x io^ (i+|V|/4.89)N e /2] 

(focused ion collection onto structure) 

J” = [1.2 x io _1 5 (i+|V|/4.89) l/2 N e /2] 

(focused ion collection onto array masts) 
and A ISS = 20 m 2 , A mast = 7.5 m 2 , and A sa = 6.3 m 2 . 

In 4 > = 5.194 - 0.02415 V[(N e /T e ) l/2 ] (2 arrays) 

In c|y = 5.2147 - 0.01741 V[(N e /T e ) l/2 ] (4 arrays) 


Fitting the Difference of Two Exponentials 

V ob s = -cf> [e at - e bt ] , because both electron and ion collection 
decay with time because of surface charging. 

tp k = (In a -In b)/(a-b). 

Values of “a” from our fits ranged from 0.43 to 5.5 s 1 , (for l/e 
times of 0.18 to 2.3 seconds) with a mean of 1.61 and median of 
1.24 s" 1 (0.62 and 0.80 seconds, respectively). 

Values of “b” from our model fits ranged from 0.056 to 0.60 s-i 
(for l/e times of 1.7 to 18 seconds), with a mean of 0.21 and 
median of 0.18 s-i (4.8 and 5.6 seconds, respectively). 

Values of the ratio of a/b ranged from 2.40 to 52.7, with a mean 
of 10.6 and median of 6.3. 

Finally, the amplitude cf> flt ranged from 7.35 to 91.9 volts, with a 
mean of 28.7 and median of 23.2. 


RCE fit by the difference between two exponentials 




1st exponential time constant 
a (s' 1 ) vs Ne (ram ion charging on structure elements) 

a(s 1 ) vs Ne 



Ne (10 11 rrv 3 ) 


M«-e5cp onential time consta 
b (s 1 ) vs V(N e T e ) (thermal electron charging on arrays) 



;i 3U5 


i = theory 


(p vs V (Ne/Te) 



V (Ne/Te) 




RafrnaTized AmplTtude^/ersOi^^^Karr^ngle 
(line = cosine law) 
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4 active arrays 


<p vs V (Ne/Te) 



♦ Phi (fit) 

— Phi' 4 arrays 
— Phi 2 arrays 
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